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ABSTRACT 



We present infrared observations of the ultra-compact H II region W3(OH) made by 
the FORCAST instrument aboard SOFIA and by Spitzer/IRAC. We contribute new 
wavelength data to the spectral energy distribution, which constrains the optical depth, 
grain size distribution, and temperature gradient of the dusty shell surrounding the 
H II region. We model the dust component as a spherical shell containing an inner 
cavity with radius ~ 600 AU, irradiated by a central star of type 09 and temperature 
~ 31,000 K. The total luminosity of this system is 7.1 x 10 4 L©. An observed excess 
of 2.2 — 4.5 /im emission in the SED can be explained by our viewing a cavity opening 
or dumpiness in the shell structure whereby radiation from the warm interior of the 
shell can escape. We claim to detect the nearby water maser source W3 (H2O) at 31.4 
and 37.1 /jm using beam deconvolution of the FORCAST images. We constrain the 
flux densities of this object at 19.7 — 37.1 /im. Additionally, we present in situ obser- 
vations of four young stellar and protostellar objects in the SOFIA field, presumably 
associated with the W3 mo lecular cloud. Results from the model SED fitting tool of 



Robitaille et alJ (120061 : 120071 ) suggest that two objects (2MASS J02270352+6152357 and 



2MASS J02270824+6152281) are intermediate-luminosity (~ 236-432 L Q ) protostars; 
one object (2MASS J02270887+6152344) is either a high-mass protostar with luminosity 
3 x 10 3 L or a less massive young star with a substantial circumstellar disk but depleted 
envelope; and one object (2MASS J02270743+6152281) is an intermediate-luminosity 
(~ 768 Lq) protostar nearing the end of its envelope accretion phase or a young star 
surrounded by a circumstellar disk with no appreciable circumstellar envelope. 
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1. INTRODUCTION 



Ultra-compact H II regions can be found throughout the galaxy surrounding newly formed 
massive stars just completing the accretion stage and entering their main sequence lifeti mes. UCHs 
are s mall (< 0.4 pc), hot (> 100 K), and massive (> 10 3 M Q ), and emit 10 4 - 10 6 L (jChurchwell 
20020 . These regions of ionized gas are surro unded by molecula r and dust clouds, out to as far as 
ten times the radius of the UCHs themselves (jConti et al.l 120071 ). This causes attenuation of much 
of the emitted luminosity, as it is absorbed and reradiated in the infrared. 

W3(OH) is one of the largest and best studied UCH II regions known in the Galaxy. It is at a 
distance of approximately 2.04 kpc in the Perseus arm (?). Like most UCH II regions, is surrounded 
by its natal dust and molecular gas envelope as well as a larger giant molecular cloud encompassing W3(OH) 
and W3(H 2 0), W3 Main, and AFGL 333; these are all regions of potentially triggered sta r formation based 



on their positions on the dense outskirts of a much less dense cavity in the W3 G MC (Ruch et all 



Moore et al]|2007l ). Further review of star formation in the W3 GMC is presented bv lMegeath et al.l ([2008) 



2007 



Line emission from the molecular gas surrou nding W3(OH) has been mapped at sub millimeter and ra 



dio w avelengt hs in the transitions of OH and H9O (IMader. Johnston & Moran 

1984). NHq Jwilson. Bieging fc Downesll978l:IZeng et alll984tlTieftrunk et al.lll998f ). CIUOH (jMenten et al 
19921 ). C 18 Q jwink et al.lll994h . and C 17 Q jWvrowski et al.lll997h . This molecular emission is found out to 



1978), HCN dTurner fe W elch 



a 1' diameter region. These studies have particularly focused on further characterizing the molecular gas 
and the OH and methanol masers within and around W3(OH), as well as the hot core/ water maser region 
W3(H 2 0) approximately 6" to the east. In continuum emission, however, the dust cocoon surrounding the 
UCH II region comes into primary focus. 



The dust component to the UCH II W3(OH) was detected by I Wvnn- Williams. Becklin. fc Neugebauer 

(J1972J) at ~ 1 — 20 /im; the dust component is optically thick in the near- and mid-infrared. ? studied the 



cold dust at wavelengths of 350 /im and 1.3 mm. Using a spherically symmetric radiative transfer model, 
they concluded that the inner cavity (< 2 x 10 17 cm) of the UCH II region is depleted of dust, rather 
than having dust density that increases approaching the central star. A large amount of visual extinction 
(A v ~ 67) in a thick outer shell wa s required to explain the decline in emission at A < 5 /im. Utilizing 
airborne data, ICampbell et al.1 (|l989n showed that the dust cocoon is optically thick in the far-infrared and 
developed a more detailed model containing an H II region and a cavity. Surrounding the cavity is a dusty 
region with a free-fall density distribution and a temperature gradient through the dust cloud. More recently, 
Stecklum et al.l (|2002f) presented high spatial resolution, ground-based 10 and 20 /im images of W3(OH). 
From these data, the model for the dust shell was further developed and contained a Gaussian density 
distribution with an inner cavity of radius 2270 AU and stellar luminosity of 8 x 10 Lq. 

In this work, we present new, high spatial resolution observations of t he W3(OH) region in the wave- 
length range ~ 3. 6 - 40 fim obtained with the FORCAST instrument faerter et al.l l2012h on SOFIA 
( Young et al 2IJ12}) and with the Infrared Array Camera (IRAC; Fazio et al. I2OO4J ) on the Spitzer Space 
Telescope (j Werner et al.l 120041 ) . These wavelengths are critical for determining the spectral energy distri- 
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bution (SED) of the W3(OH) d ust cocoon and thus fo r making a measurement of its total luminosity. We 



combine our data with 2MASS jSkrutskie et al.ll2006h fluxes and other data published in the literature to 



construct SEDs for the W3(OH) dust component. We model the dust component as a dusty shell around 
the H II region and compute the emergent flux using radiative transfer code. In addition, we present in 
situ observations of four stellar and protostellar objects in the SOFIA/FORCAST field. We fit the SEDs 
of these objects to those of high- and intermediate-mass protostars containing dusty circumstellar envelopes 
and circumstellar disks. We discuss the physical properties of all these objects. 



2. OBSERVATIONS 



2.1. SOFIA/FORCAST Observations 

Observations were conducted on December 8, 2010 aboard the Stratospheric Observatory for Infrared 
Astrono my (SOFIA), using the Cornell-built Faint Object Infrared Camera for the SOFIA Telescope (FOR- 



CAST) (|Adams et al 



120101 . FORCAST uses a Si:As BIB detector array at wavelengths A < 26//m and a 
Si:Sb BIB array at A > 26/im. The detector has a 0.7638" pixel scale (rectified) over a total field of view 
of 3.4' x 3.2'. Four filters were utilized for these observations, with central wavelengths (and bandpasses) 
at 19.7 (5.5), 24.2 (2.9), 31.4 (5.7), and 37.1 (3.3) fim. The dichroic beamsplitter was used to obtain the 
wavelength pairs 19.7/37.1 and 24.2/31.4 /j,m. Asymmetrical chopping with respect to the optical axis was 
applied in order to avoid telescope coma in the on-source beam. The c hop throw was 5' an d was designed to 
chop off nearby nebulosity. A 5-point dither pattern in C2NC2 mode (|Herter et al.ll2012l ) was implemented 
to remove bad pixels during post-processing. Integration times were 30 sec at every dither position which 
yielded a total of 150 se c in each filte r confi guration. The data were pipeline processed with the reduction 
algorithms described in iHerter et al.l (|2012 ). A global solut ion to the calibration from the SOFIA Early 
Science phase was applied to the images (jHerter et al.ll2012l ). We discuss flux extraction of sources in Sj3] 
Finally, we a pplied color correct ions to the extracted fluxes as computed by an instrument model with model 
atmosphere (jHerter et al.ll2012 ). 



2.2. Spitzer/IRAC Observations 



We used data obtained with IRAC on Spitzer. The data were taken from observations that were 
obtained in high dynamic range (HDR) mode, whereby two images are taken in succession at 0.4s and 
10.4s integration times. The brightest objects in the field, including W3(OH) were saturated in the longer 
frames, so we used the 0.4s exposure time frames to construct the mosaic that was used for the bright 
source photometry. The following AORIDs were used: 5050624, 19305728, 20590592, 38744064, 38757632, 
38763776, 38769408, 38770944, 38790656, and 38801408. We utilized the Basic Calibrated Data (BCD) 
version S18.5 products from the Spitzer Science Center (SSC) standard data pipeline. For the 3.6 and 4.5 
/im bands, data from the cryogenic and warm mission were combined. For the 5.8 and 8.0 /im BCDs, bright 
source artifacts such as "banding" (jHora et al.l 120041 ) we re removed using the IMCLEAN im age processing 
routine^. The images were mosaiced wi th IRACproc (Schuster. Marengo, fc Patten 2006 ) which uses a 
version of the mopex mosaicing software (jMakovoz fc Khanl 120051) developed at the SSC. The final mosaic 



http: / /irsa.ipac. caltech.edu /data/SPITZER / docs /dataanalysistools/tools / contributed/irac/imclean / 
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was made with a pixel scale of 0.863"/pixel. 



3. RESULTS 



3.1. Imagery and SOFIA/FORCAST detections 

Reduced, multi- wavelength IRAC and FORCAST images are shown in Figs. Q] and [5J The measured 
image quality in the FORCAS T images was consi stent with the overall image quality achieved by SOFIA 
during the Early Science period (jHerter et al.ll2012l ). W3(OH) is spatially resolved in the FORCAST images. 
Within W3(OH), the dual peaks in flux correspond to a bright cometary South western H II region and a 
dimmer, more eliptically-shaped Northeastern H II region (jStecklum et al]|2002n . 



We detect four mid-infrared point sources in the SOFIA field. These sources are designated in the 
2MASS catalog as J02270352+6152357, J02270743+6152281, J02270824+ 6152281, and J0227088 7+6152344. 
There are also 2 fim counterparts to these sources in the K' band image of Tieftrunk et al. ( 1998^ . However, 
there is a set of rich clusters of low mass stars in this region, so confusion with neighboring sources precludes 
our making an accurate measurement of their 2 fim fluxes. In Fig. [3J we label the sources with their 
corresponding designations. All four sources were detected by SOFIA/FORCAST at 19.7, 24.2, and 31.5 
/im; J02270352+6152357, J02270824+6152281, and J02270887+6152344 were also detected at 37.1 fim. 
In addition, all four sources were detected in the Spitzer images at wavelengths of 3.6, 4.5, 5.8, and 8.0 
fim; however, the 8.0 fim detection for J02270743+6152281 suffers from contamination by relatively bright 
nebulosity extending from below the source. The fluxes for J02270743+6152281, J02270824+6152281, and 
J02270887+6152344 were extracted from FORCAST images by fitting Gaussian functions to flux line profiles 
and from Spitzer images by aperture photometry. The corresponding flux densities are given in Table [TJ 
Flux extraction for J02270352+6152357 is discussed in 



3.2. Deconvolved FORCAST images 



We performed beam deconvolution on each of the FORCAST images in order to search for the hot 
core W3(H20) an d to measure the size of W3 (OH). The images were deconvolved using the maximum 
likelihood method ([Richardson! 1 1972t lLucvlll974j ). Like all deconvolution methods, knowledge of the point- 
spread function (PSF) of an unresolved source is needed at each wavelength. The delivered PSF can change 
due to different wind loads on the secondary mirror and differences in the telescope flexure as a function 
of telescope position. To mitigate these effects on ground-based telescopes, high S/N observations of mid- 
infrared bright stars are usually taken immediately before and/or after each science target observation and 
as close to the science target as possible (< 1° away) so as to get the best PSF calibration for use in the 
deconvolution procedure. However, in our case, finding a PSF star that is bright enough at wavelengths out 
to 40 ^tm that fulfills these requirements is nearly impossible and thus was not attempted. Standard star 
images taken throughout our flight similarly did not have sufficient S/N for use as a PSF calibrator at the 
longer wavelengths. Therefore, we used these standard stars observations to determine an average FWHM 
for each wavelength for the flight. Then artificially generated PSFs (using an Airy pattern calculated from 
the wavelength, telescope diameter, and central obscuration diameter) were constructed and convolved with 
a Gaussian to achieve PSFs with FWHMs that equaled the measured average FWHMs of the standard stars. 
These idealized PSFs were then used in the deconvolution procedure. These deconvolved images compare 
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favorably to simple unsharp masking of the original images, and hence all of the substructures revealed in 
the deconvolved images are believed to be real with high confidence. Final image resolutions are about a 
factor of two better than the natural resolutions for each image. 

We claim to detect unresolved emission from the hot core W3(H 2 0) at 31.4 and 37.1 /im. Figs. [3] 
and [TJ show the region near W3(H 2 0) in the deconvolved images at 31.4 and 37.1 /im, respectively. In each 
case, the central core of W3(OH) has been fit to a 2-dimensional Gaussian fu nction and subtracted from th e 
respective image. Figs. [3]and|3]also contain 8.4 GHz VLA observations from ( Wilner. Reid fc Menten 1999h . 



The radio continuum image shows a cometary component extending in the North east direction from th e peak 



of W3(OH). The radio image also shows emission from the Northeast H II region (jStecklum et al.ll2002l) above 



the c ometary component. In the 8.4 GHz continuum, the W3(H20) clumps A, B, and C (jWyrowski et al 



1999) are resolved. Figs. [3] and 0] show residual emission at 31.4 and 37.1 /im in the vicinity of W3(H20). 
Moreover, the position of this residual emission shifts in wavelength from close to W3(H20) C at 31.4 jam 
towards W3(H 2 0) A at 37.1 /im. 

Flux extraction for W3(OH), W3(H 2 0), the Northeast H II region, and 2MASS J02270352+6152357 
was performed using the deconvolved images. The resolved diameter of W3(OH) in the deconvolved images 
(Fig. [5]) is approximately 40 pixels, corresponding to ~ 63000 AU. We used a 15-pixel aperture radius to 
extract the flux density at 19.7 /im and 20-pixel aperture radii to extract the flux densities at 24.2, 31.4, 
and 37.1 /im; we then subtracted the flux densities measured for the Northeastern H II region and 2MASS 
J02270352+6152357 (discussed next). The resulting integrated flux densities for W3(OH) are listed in Table 

HI 

Fig. [SJshows boxed regions where the flux profiles for W3(OH) and W3(H 2 0) were extracted (integrated 
vertically) and fit to Gaussian functions (Fig. [6]). Each of the profiles shows a secondary peak coincident 
with the location of W3(H 2 0). The area under the profile fit to W3(H 2 0) yields its flux. The extracted flux 
densities for W3(H 2 0) are also listed in Table [TJ We discuss our detection of this source further in §4.21 
The flux densities for W3(OH) measured in this fashion agree to within ~ 10% of the flux density measured 
from large aperture photometry, which is listed in Table [TJ 

Fig. [5] also shows the regions selected for line profiles of the Northeastern H II region (integrated 
horizontally) and 2MASS J02270352+6152357 (integrated vertically). The latter source is flagged in the 
2MASS catalog as confused with neighboring objects; thus we do not consider its 2 /im flux. The line 
profiles for these 1-dimensional integrations are shown in Fig. [6] and their flux densities, determined by 
Gaussian line fits, are listed in Table [TJ 



DISCUSSION 



4.1. W3 (OH) 

In Fig. [7| , we show th e SED for W3 (OH) with additional data taken from t he lite rature: 2MASS (2.2 



urn UStecklum et all (120021) (8.8. 12.7. and 17.8 /im), MSX (21 /im, Egan et al. 119991 ). IRAS (60 and 100 



/im), and Chini. Kriigel fc Krevsa ( 19861 ) (1.3 mm). The IRAS flux densities were taken as upper limits due 
to contamination from nearby sources. 



We model the dust component as an optically thick, dusty shell around the H II region, irradiated at 
the inner boundary by the central star with surface temperature 3.11 x 10 4 K. The emergent SED of the 
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model was computed using the DUSTY radiative transfer c odjl. The densit y distribution of this model 
(p oc r~ p ; p = 1.5) is based on a free-fall density profile (e .g. Hartmann 20091) . We chose the composition 



of the dust to be 53% silicates and 47% graphite grains ( Draine fc Lee "Il984 ) . We consider a grain size 
distribution n(a) oc a -9 , whereby q = 3.5 (Mathis, Rumpl, & Nordsieck [1977J; hereafter MRN), and with a 



minimum grain size of 0.001 /jm and a maximum grain size of 0.25 /im (MRN; Sellgren ll984f) . We set the 
temperature of the inner edge of the shell at 400 K for this model. This is substantially cooler than the 
dust sublimation temperature (~ 1500 K), indicating the presence of a large cav ity depleted of dust and 
consistent with previous work ( Chini. Kriigel fc Krevsa 19861 : Stecklum et al. 2002 ). The total luminosity of 
this model is ~ 7.1 x 10 4 Lq. The SED of this model is shown in Fig. [7] as the dotted line and the relevant 
model parameters are listed in Table [5J 

Although they both use an optically thick, free- fall de nsity shell in the IR, this model differs quanti- 
tatively from the model presented in Stecklum et al. ([20021 ) . The central star is of later O-type and cooler 
surface temperature. The inner dust shell radius is nearly four times smaller and the outer radius nearly 
twice as compact as in Stecklum et al. These parameters are necessary in order to explain the mid-IR 
emission in the range 17 — 37 /im. This work finds that the total luminosity is slightly lower than previous 
estimates. 

Note this model cannot account for the excess emission at 2.2 — 5.8 /im. One explanation for this 
emission is dumpiness in the dust cloud, or a cavity opening, which allows short wavelength radiation from 
warm dust to leak through holes in the cloud. We model radiation from the warmer interior of the shell 
as a single-temperature (425 K) blackbody component. The plausibility of this model is evident in Fig. [71 
whereby an excellent fit to the data is aChieved with a temperature consistent with that of the interior of 
the dust shell. There may also be a scattered light component at these wavelengths, but we do not model 
such a component since its geometry is relatively unconstrained. 

An alternate explanation for the excess at 2.2 — 5.8 /im is a relatively high abundance of very small 
grains compared with large grains. In order to fit the data at 2.2 — 5.8 /im, we would need to consider a 
model with a modified MRN grain size distribution (q = 4.7). However, in the ionized region, small grains 
can become super-heated and destroyed. The emission would come from a relatively high abundance of 
small grains in the infalling shell; such grains would need to be primordial. Given the jets that are seen 
in molecular emission, it is most likely that the excess short wavelength radiation originates in a cavity 
opening. The prediction one can make from this explanation is the presence of continuum emission in the 
amporphous silicate absorption band at 9.7 /Ltm, which decreases the depth of the absorption feature. Mid-IR 
spectroscopy in the 8 — 13 /im window could be used to test this prediction. 



4.2. W3(H a O) 



In the radio continuum, the hot core W3(H2Q ) consists of 3 clumps designated as A, B, and C 
(jWvrowski et alJll99a Iwilner Reid fc Menta3ll999l ). An early claim to a detection of W3(H 2 0) in the 
infrared came from iKeto et al.l (|l992T) who presented ground-based observations at 12.2 /im showing emis- 
sion at the location of W3(H2 0) C, which lies approximately 4" East of W3(OH). No emission was seen from 
sources A and B, which lie approximately 6" East of W3(OH). However, this observation was not substan- 



7 Ivezic, Z., Nenkova, M., and Elitzur, M. 1996, University of Kentucky, 



http://www.pa.uky.' 
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tiated by IStecklum et al.l (|2002l) who did not detect any emission from these sources at 8 — 12 /Ltm, despite 
their spatially resolving W3(OH). In the FORCAST data, the separation of the peaks in the double-peaked 
line profile (Fig. EJ yields the separation between W3(H 2 0) and W3(OH), which we find to be 5.0" at 37. 1 
/im and 4.4" at 31.4 /im. Components A an d C form a massive p rotobinary system (jMinh fc Chenl 12007). 
Based on the gas chemistry of these sources, iMinh fc Chenl (|2007l ) suggest that component A may be more 
deeply embedded and younger than component C within < 10 4 yr. This would be consistent with our results 
that show a wavelength dependence for the 19.7 — 37.1 /Ltm emission. However, it should be noted that this 
dependence on wavelength may also result from combined temperature and optical depth gradients. 



4.3. Protostars and Young Stars with Disks in the SOFIA Field 



We construct SEDs for the four point sources in the FORCAST field and co mpare them with model 



SEDs of protostellar and young stellar objects using the online SED fitting tool of Robitaille et al. ( 20061; 



20071 ). The SEDs and those of the best-fit models are shown in Fig. [8j We provide comments on each object 
in gXU- JOT] 



4.3.1. 2MASS J02270352+6152357 

The best fitting model for J02270352+6152357 corresponds to a protostellar object with a mass of 3.96 
M and luminosity of 236 L Q . The model indicates that this object is young and cool, with an age of just 
under 5000 years and a temperature in the low 4000s of Kelvins. It is undergoing active envelope accretion 
from a massive envelope (11.7 M©) onto a low-mass disk (3.92 x 10~ 3 ). The self-consistency of the top 
ten models with this best-fitting model suggest that this object is indeed a young, intermediate-luminosity 
protostar. 



4.3.2. 2MASS J02270743+6152281 

The SED fitting tool produced two families of best-fit models for J02270743+6152281, with the first 
family headed by the best-fitting SED and the second by the second-best fit. The families differ in the 
object's stage of envelope accretion and disk formation, with one family representing an object still actively 
undergoing envelope accretion and the other corresponding to a more highly-developed disk with a depleted 
envelope. 

The best fitting (total x 2 = 82.14) model, representing the first family of models within the top ten 
best fits, indicates that this object is a protostellar object with a mass of 5.69 M© and a luminosity of 768 
Li©, surrounded by a relatively substantial envelope (0.153 Mq) and a smaller disk (3.36 x 10~ 2 M©). This 
model has a low but nonzero rate of envelope accretion, and suggests an age of 6.47 x 10 5 yr, with an internal 
temperature of 16,500 K. These parameters point to an intermediate-luminosity protostar nearing the end 
of its envelope accretion phase. 

The second-best fit (total % 2 = 82.72), and correspondingly the second family of fits, represents an 
object with a fully formed disk and a completely depleted envelope no longer accreting material onto the 
disk. This model protostar has an intermediate mass somewhat higher than the alternative family, of 9.17 
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M . Its envelope is negligible (10 -5 M©), and its disk has a mass of 2.1 x 10~ 2 M Q . It suggests an age 
of just over 10 6 yr and a temperature of 24,500 K. Given the protostellar mass and these parameters, this 
family of models likely corresponds to a young star with a circumstellar disk. 

The two families of models differ primarily in the presence or absence of an amorphous silicate absorption 
feature at 9.7 /im, which is typically seen in an envelope-dominated SED. Highly sensitive mid-infrared 
spectroscopy could be used to resolve the degeneracy between these two families of models. 



4.8.8. 2MASS J02270824+6152281 

Model fitting for J02270824+6152281 presents a convergent set of parameters. The best fit model 
corresponds to a 6.11 M Q protostar with a luminosity of 432 L©. This model suggests the object is a young 
protostar (approximately 8000 years) with a temperature of around 4000 K. The object is embedded in a 
large 35.1 M© envelope with a disk of mass 1.48 x 10 -2 M©, indicating that this object is likely experiencing 
ongoing disk formation. The relatively high envelope accretion rate displayed by this model (1.39 x 10~ 3 M© 
yr -1 ) supports this conclusion. The models predict substantial far-IR emission arising from the envelope, 
meaning followup observations in the far-IR range could confirm our assessment. 



4.3.4. 2MASS J02270887+6152344 

For J02270887+6152344, the data again result in two possible families of models, representing cither 
an intermediate- to high-mass protostellar object with a large envelope or an older intermediate-mass young 
star with a clearly defined disk and minimal envelope. The top three fits all corresponded to the former of 
these models, indicating its greater likelihood for accuracy. 

The best-fitting (total \ 2 = 15.03) model for this source represents a 10.5 M© protostar surrounded 
by a natal dust envelope containing 240 M© of material and a negligible disk. The model object has a 
luminosity of 3000 L©. An age of approximately 2000 yr and a temperature in the 4000s of Kelvins indicate 
that the object is protostellar in nature. 

The second group of well-fitting models is represented by the fourth-best fitting (total x 2 = 34.34) 
model. This model corresponds to a young star of mass 9.8 M©, with a depleted envelope (mass approxi- 
mately 2 x 10~ 5 M©, indicating that the envelope has nearly entirely accreted onto the disk). This model 
possesses an age of over 10 6 years and an internal temperature of 25,400 K, indicating a more advanced stage 
in development than the first family of models presented. 

These families of models diverge at far-infrared wavelengths. With only present data, we therefore 
cannot rule out either the higher-mass, younger protostellar object or the intermediate-mass, more evolved 
young star. Further observations of this object at far-IR and/or submillimeter wavelengths are required to 
resolve the degeneracy in the model parameters. 
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5. CONCLUSIONS 

We present SOFIA/FORCAST and Spitzer/IRAC observations of the UCH II region W3(OH) in the 
wavelength range 3.6 — 37.1 /im. These data, combined with other published data, have been used to 
constrain the optical depth, grain size distribution, and temperature gradient in the dusty shell surrounding 
the H II region. The total luminosity of W3(OH) is 7.1 x 10 4 L Q , indicating that the central star is an 09 
star with surface temperature ~ 31,000 K. A clumpy dust distribution or cavity opening revealing warm 
interior grains is necessary to explain excess emission at 2.2 — 4.5 /im. 

We detect the hot core W3(H20) at 31.4 and 37.1 /zm, and constrain its flux density at 19.7 — 37.1 /Ltm 
using deconvolved FORCAST images. 

In addition, SEDs have been constructed for four you ng stellar or protostel lar objects which lie in 
the SOFIA/FORCAST field. The model SED fitting tool of iRobitaille et all |2006l ) was used to determine 



the nature of these objects. 2MASS J02270352+6152357 is an intermediate-luminosity protostar undergoing 
envelope accretion; 2MASS J02270824+6152281 is most likely a very young intermediate-mass protostar with 
a large natal envelope; 2MASS J02270887+6152344 is a high-luminosity object which is either a protostar 
with ongoing envelope accretion onto a young disk or a young star with a circumstellar disk and a depleted 
envelope; and 2MASS J02270743+6152281 could be an intermediate-luminosity protostar or potentially a 
young star with a developed disk and an almost entirely depleted envelope. Further observations in the mid- 
IR, far-IR and/or submillimeter range(s) are required to definitively characterize 2MASS J02270887+6152344 
and 2MASS J02270743+6152281. 
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Fig. 1. — False-color IRAC image of the W3(OH) region at 3.6 /xm (blue), 5.8 /j,m (green), and 8.0 
/xm (red). The box indicates the region of the SOFIA field that is shown in Fig. [2l 
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Fig. 2. — False-color image of the W3(OH) region at 3.6 /im (blue, IRAC), 19.7 //m (green, FOR- 
CAST), and 37.1 /im (red, FORCAST). The image at 3.6 /im was convolved with a Gaussian kernel 
to match the spatial resolution of FORCAST at 37.1 /um. The positions of W3(OH) and the four 
point sources that were detected by SOFIA/FORCAST are indicated. 
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Fig. 3.— Deconvolved SOFIA/FORCAST image of W3(OH) at 31.4 /xm, scaled linearly in flux. 
The image is centered on the coordinates RA=02 h 27 m 03 s .87 and DEC=+61°52'24".6 (J2000). The 
core of W3(OH) has been subtracted using a 2-dimensional G aussian fit to its radia l brightness 
profile. The locations of sources A, B, and C in W3(H20) fro m I Wvrowski et al.l (119991 ) are s hown. 
The green contours represent 8.4 GHz VLA observations from lWilner. Reid Mentenl (|1999i ). The 
contours levels are 0.00004, 0.0001, 0.00018, 0.0003, 0.0005, 0.001, 0.005, 0.01, 0.015, 0.016, 0.02, 
0.022, and 0.03 Jy/beam. An assumption is that the peak of the 8.4 GHz contours coincides with 
the peak of the 19.7 — 37.1 /xm emission. Also seen is the H II region to the Northeast of the peak 
of W3(OH). 
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Fig. 4. — Same as Fig. [3] for 37.1 /im. 
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Fig. 5.— SOFIA/FORCAST deconvolved images at 19.7, 24.2, 31.4, and 37.1 fim, scaled logarith- 
mically in flux with rain bow color intensity m apping. The locations of W3(OH), W3(H20), the 
Northeastern H II region (jstecklum et al.ll2002l ). and the point source 2MASS J02270352+6152357 
are indicated in the upper left panel. The boxed regions were the regions selected for flux extrac- 
tion; 1-dimensional integrations were performed over the short dimensions of the boxes and the 
resulting line profiles are shown in Fig. [6j 
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Fig. 6. — 1-dimensionally integrated and normalized line profiles (points) and best Gaussian fits for 
W3(OH) (magenta line), W3(H20) (blue line), the Northeastern H II region (green line), and the 
point source 2MASS J02270352+6152357 (red line). The 1-dimensional integrations were performed 
along the short directions of the boxes in Fig. [5j The black lines represent the sum of the fits. The 
extracted flux densities for all sources are given in Table HJ 
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Fig. 7. — SED of W3(0H) and model fits. Filled circles: Combined data points (see text for 
references). The IRAS fluxes were taken as upper limits, indicated as downward arrows at ~ 60 and 
100 /an. Dotted line: DUSTY model SED. Dashed line: Improvised single-temperature blackbody 
originating from interior, warm dust emission through a cavity opening. Solid line: Total SED 
from DUSTY model and improvised single-temperature blackbody. 
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Fig. 8. — SEDs of the four point sources in the SOFIA field (fill ed circles) and the be s t-fit p roto- 
stellar SED model (solid line) from the online SED fitting tool of iRobitaille et all fcood : 120071 ) . For 
2MASS J02270743+6152281 and 2MASS J02270887+6152344, the best-fitting alternative, disk- 
dominated model SEDs are also shown (dashed lines). 



Table 1. Positions and flux densities in Janskies for W3(OH), W3(H20), and four mid-IR sources in the SOFIA field. 
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Table 2. DUSTY parameters for best-fit model of the W3(OH) dust component. The 
parameters are stellar temperature inner shell radius Ri n ; outer shell radius R ou t, inner shell 
boundary temperature Tj n , outer shell boundary temperature T out , grain composition, density 
distribution parameter p, grain size distribution parameter q, and optical depth T37 at 37 fxm.. 



Parameter 


Value 


T* 


31145 K 


Rin 


576 AU 


Rout 


29942 AU 


J- in 


400 K 


Tout 


26 K 


Composition fraction, silicates 


53% 


Composition fraction, graphite 


47% 


P 


1.5 


q 


3.5 


T37 


2.8 



